Control the entanglement of two atoms in an optical 

cavity via white noise 



Jing-Bo Xu and Shang-Bin Li 
Zhejiang Institute of Modern Physics and Department of Physics, 
Zhejiang University, Hangzhou 310027, People's Republic of China* 



Abstract 

Two two-level atoms within a leaky optical cavity is driven by two in- 
dependent external optical white noise fields. We investigate how entan- 
glement between two atoms arises in such a situation. The steady state 
entanglement of two atoms is also investigated. A stochastic-resonance-like 
behavior of entanglement is revealed. Finally, the Bell violation between 
atoms is discussed. 
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I. INTRODUCTION 

Quantum entanglement plays a crucial role in quantum information and quantum com- 
putation [1]. Entanglement can exhibit the nature of a nonlocal correlation between 
quantum systems that have no classical interpretation. However, real quantum systems 
will unavoidably be influenced by surrounding environments. The interaction between 
the environment and quantum systems of interest can lead to decoherence. It is there- 
fore of great importance to prevent or minimize the influence of environmental noise 
in the practical realization of quantum information processing. In order to prevent the 
effect of decoherence, several approaches have been proposed such as quantum error- 
correcting approach [2] or quantum error-avoiding approach [3,4]. 

Instead of attempting to shield the system from the environmental noise, Plenio 
and Huelge [5] use white noise to play a constructive role and generate the controllable 
entanglement by incoherent sources. They showed that the noise-assisted entanglement 
exhibits the stochastic resonance behavior. Similar work on this aspect has also been 
considered by other authors [6-8] . In this paper, we study the quantum system in which 
two two-level atoms within a leaky optical cavity is driven by two independent external 
optical white noise fields. We investigate how entanglement between two atoms arises 
in such a situation. It is shown that white noise exhibits dual aspects, i.e., playing ei- 
ther a destructive or a constructive role in quantum information processing. Recently, 
Clark and Parkins [9] proposed a scheme to controllably entangle the internal states 
of two atoms trapped in a high-finesse optical cavity by employing quantum-reservoir 
engineering. By making use of laser and cavity fields to drive two separate Raman 
transitions between stable atomic ground states, the two atoms is effectively coupled 
to a squeezed reservoir. Phase- sensitive reservoir correlations leads to entanglement 
between the atoms. Different from their scheme, we will focus here on the problem 
of generating entanglement when only incoherent sources are available and show that 
controllable entanglement can arise in this situation. We show that, if two atoms are 
simultaneous driven by two independent white noise field with the same intensity, the 
entanglement between them is suppressed and eventually completely destroyed by the 
noise. However, in another case, in which only one atom is exposed in white noise field, 
the steady state entanglement of the two atoms is non-monotonic function of both the 
intensity of noise driving field and the spontaneous decay rate. A double resonance 
behavior emerges. Moreover, the threshold value of the spontaneous decay rate, below 
which there is not any steady state entanglement, increases with the intensity of noise 
field. 

This paper is organized as follows: In section II, we outline the experimental set 
up, in which two atoms are trapped in a optical cavity and driven by the thermal field. 
The spectral width of the thermal field is large compared to the linewidth of the atomic 
transition so that its effect is that of a white noise source, we model this system by 
a master equation and give a explicit analytical solution of the time evolution density 
matrix. In section III, based on the density matrix, we obtain the analytical expres- 
sion of the concurrence characterizing the entanglement between two atoms. Both the 
entanglement during the time evolution and the steady state entanglement are inves- 
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tigated. A conclusion is given in section IV. 



II. THE MASTER EQUATION DESCRIBING TWO 
ATOMS TRAPPED IN A OPTICAL CAVITY AND 
DRIVEN BY NOISE FIELD 

The system we consider here is two atom trapped in a optical cavity. The atoms are 
driven by two independent thermal fields and separated by a large enough distance 
that they feel no direct dipole-dipole interaction. The cavity has a field decay k and a 
frequency uj. The two eigenstates of the individual atom (|0), |1)) constitute the qubit 
states. The master equation for the total system density operator is (h = 1) 

^ = ~i[H, p] + CcavP + jC at p, (1) 

where 

2 2 

H = ut Ja + f - |0)a-<0|) + siyiOU 1 ! + a|%(0|), (2) 

Z 3=1 3=1 

where a and a* are the annihilation and creation operators of the cavity field with 
frequency uj, and uj is the transition frequency of the atoms and g is the atom-cavity 
coupling constant. The Liouvilleans C cav p and C at p are given by 

£cavP — K>(2apa — wap — pa* a), (3) 

and 

CatP = E(^t + i)r^2|oui| P |iuo| - |i)„(i| P - p|i)^(i|) 

3=1 

+ E^ ) rW(2|l) j ,(0|p|0),,(l| - |0^(0|p-p|0),,(0|) (4) 

3=1 

where describes the coupling strength of the jth atom to the external fields and 
n 0')p(i) j s transition rate due to the thermal field. The spectral width of the ther- 
mal field is large compared to the linewidth of the atomic transition so that its effect 
is that of a white noise source. Here, rij) can be interpreted as an effective photon 
number and that spontaneous decay of the atom out of the cavities is included in this 
scenario via the + 1 term. 

In the large detuning limit, i.e., A = uj q — uj ^> gy/n + 1 with n being the mean 
photon number of the cavity field, there is no energy exchange between the atomic 
system and the cavity. We can obtain the effective Hamiltonian H e [10,11] 

2 2 

He = ^[£(11)^(11^ - |0)#<0|ata) + |l)n(0| ® |0) 22 (1| + |0>u<l| ® |1> 22 (0|]. (5) 
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+ ^p s (|l>n(l| + |1) 22 (1| + |l)n(0| ® |0)22(1| + |0) n (l| ® |1) 22 



+ E("£ ) + i)r^(2|o)„(i| Ps |i)„(o| - \i) Sj (i\ Pa - Ps\i)M 



The first and second terms describe the photon number dependent Stark shifts, and the 
third term describes the dipole coupling between the first and second atoms induced by 
the virtual photon process. When the cavity mode is initially in the vacuum state |0), 
it will remain in the vacuum state throughout the procedure, the effective Hamiltonian 
reduces to 

H e = ^(|l)n(l| + |1) 22 (1| + |l)n(0| <g> 10)22(11 + |0) n (l| ® |1)22(0|). (6) 

As the cavity mode will then never be populated, we can disregard the cavity mode in 
the following. Now, the master equation (1) can be reduced to 

^ = -i^(|l) n (l| + |1) 22 (1| + |l) n (0| <g> |0) 22 (1| + |0) n (l| <g> |l) 22 (0|)p s 
f 

'A' 
2 

3=1 

+ ^n^)(2|l) J ,-(0|p s |0),,(l|-|0)^(0|p s - Ps |0),,(0|), (7) 
j'=i 

where p s is the density matrix describing the subsystem containing only the two atoms. 
Firstly, we discuss the case with n>P = = n T and = = V, i.e., two atoms 
are driven by two independent thermal fields with the same intensity. We assume that 
the atom 1 and atom 2 are initially in the pure product state |l)i <E> |0) 2 . Then, the 
explicit analytical solution of the master equation (7) can be obtained as follows, 

p s (t) = pn(*)|l)n(l| ® |1) 22 (1| + p 22 (0|l)n(l| ® |0) 22 (0| 

+P3 3 (0|0)ll(0| |1)22(1| +P44(*)|0)ll(0| ® |0) 22 (0| 

+p 23 (0|l)n(0| ® |0)22(1| +P3 2 (t)|0)n(l| ® |1) 22 (0| (8) 

Pn(t) = f0 W ^, 2 K + e-( 4 - +2 ) r * - (n T + l)e~^ + ^% 
{2n T + 1) 

,2 11 0„2j 



with 



p 41 (t) = " r +|, 2 [ nT + 1 - e-<*-+ 2 >" - n T e-^*V'}; 
{2n T + 1J 2 

p 23 (t) = V( 4 - +2 ) r 'sin^, 

p 32 (t) = - l -e-^^sm 2 -^. (9) 



In order to quantify the degree of entanglement, we adopt the concurrence C defined 
by Wooters [12]. The concurrence varies from C = for an unentangled state to 
C = 1 for a maximally entangled state. For two qubits, in the "Standard" eigenbasis: 
|1) = 1 11) , |2) = |10), |3) = |01), |4) = |00), the concurrence may be calculated 
explicitly from the following: 

C = max{Ai - A 2 - A 3 - A 4 , 0}, (10) 

where the Aj(i = 1, 2, 3, 4) are the square roots of the eigenvalues in decreasing order of 
magnitude of the " spin- flipped" density matrix operator R = p s (o- y <g> a v )p*{a v <8> a y ), 
where the asterisk indicates complex conjugation. 

It is straightforward to compute analytically the concurrence for the density matrix 
p s (t), and the concurrence C s (t) related to the density matrix p s (t) is obtained as 
follows 

C 8 (t) = 2max{0, \p 23 (t)\ - y/ pn{t) p u (t)} , (11) 

where \x\ gives the absolute value of x. In Fig.l, we have plotted the concurrence C s (t) 
as a function of the time t and the intensity of the thermal field n T with g 2 / A = 0.2 
and T = 0.01. It is shown that the entanglement between two atoms decreases with 
%, and there is not any entanglement arising between two atoms during the time 
evolution when ut is beyond a threshold value depending on the coupling constant T 
and effective Rabi frequency g 2 / A. In Fig. 2, the concurrence C 3 (t) is plotted as the 
function of the time t and the coupling constant T of the atoms and the external field 
(Note that T is equivalent to the spontaneous emission rate if ut = 0) with two different 
values of effective photon number ut of the thermal field. In the case with ut = and 
g 2 / A = 0.2 (see Fig.2(a)), the entanglement of two atoms always arises during the time 
evolution even in the presence of atomic spontaneous emission. However, in the case 
with n T = 0.3 and g 2 / A = 0.2 (see Fig. 2(b)), a threshold value of T is found, beyond 
which there is not any entanglement arising during time evolution. All of the above 
discussions indicate that the two equal intensity independent thermal fields suppress 
the entanglement generation. But, that is not the full aspects concerning the role of 
thermal field played in the entanglement of two atoms. In the following section, we 
consider the situation in which, only one of the atoms is driven by the external thermal 
field. A different aspect of the thermal field will be found. 

III. THE STEADY STATE ENTANGLEMENT OF 
TWO ATOMS 

In the above section, we have discussed the case, in which the two atoms are simulta- 
neously driven by an external thermal field, there is not any steady state entanglement 
between two atoms in that situation. In this section, we consider the situation in which, 
only one of the atoms is driven by the external thermal field. The master equation is 
given by 

^ = -i^(|l)n(l| + |1) 2 2(1| + |l)n(0| ® |0) 22 (1| + |0) n (l| <g> |l) 22 (0|) Ps 
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q 2 

+ ^p s (|l)n(l| + |1) 22 (1| + |l) n (0| <g) |0) 22 (1| + |0) n (l| <g) |1) 22 

+ (n T + l)r(2|0) n (l|p s |l)n(0| - |l)n(l|p a - p s |l>n<l|) 

+n T r(2|l) 11 (0|p s |0) 11 (l| - |0) n (0|p s - p s |0) n (0|) 

+?7(2|0) 22 (l|p s |l) 22 (0| - |l) 22 (l|p s - p s |l) 22 (l|), (12) 

where 77 is the spontaneous emission rate of the atom 2. We assume that the two 
atoms are initially in the ground state |0)i ® |0) 2 . The explicit analytical solution of 
the steady state of the master equation (12) can be obtained as follows, 

Pst = ® |1)22(1| + P2 2 |l)ll(l| ® |0) 22 (0| 

+p s 33 \0) u {0\ ® |1) 22 (1| + P4 4 |0)ii(0| ® 10)22(01 

+ P y i)n(o| ® |o) 22 (i| + p^ 2 |o)n(i| ® 11)22(01 (13) 



with 



n 2 r 2 n 2 



P22 



P33 



(r + 77 + 2n T r) 2 (fi 2 + Tr] + 2n T Tr)) ' 

n T r[r/(r + 7] + 2n T r) 2 + fi 2 (r + 77 + n T r)] 
(r + r) + 2n T T) 2 {Q 2 + Tr] + 2n T Tr]) 

n 2 n T T(T + r] + n T T) 



P44 



(r + r] + 2n T T) 2 (n 2 + Tr] + 2n T Tr]) 

Tr](l + n T )(T + r] + 2n T T) 2 + Q 2 (T + 77 + n T lf 
(r + 77 + 2n T r) 2 (fi 2 + Tt] + 2n T r7/) 



P23 



(r + 77 + 2n T r)(fi 2 + Tr] + 2n T Tr]) ' 



^ 32 (r + r 1 + 2n T T)(Q 2 + Tri + 2n T Tr 1 y v ; 

where f2 = (? 2 /A. The concurrence C st related to the steady state p st is obtained as 
follows 



C st = 2max{0,|p 2 3]-y'p? 1 p£ 4 }. (15) 

In Fig. 3, we have plotted the concurrence C st as a function of the spontaneous emission 
rate 77 and the intensity of the thermal field nr with g 2 /A = 0.2 and T = 0.1. It is shown 
that the steady state entanglement exhibit a double stochastic-resonance-like behavior, 
which is similar with the results in Ref. [5]. The double stochastic-resonance-like be- 
havior also emerges in Fig.4, in which C s t is depicted as a function of the spontaneous 
emission rate 77 and the coupling constant T with g 2 /A = 0.2 and tit = 2. In Fig. 5, we 
plot the C s t as a function of Q and the intensity of the thermal field ny with T = 0.1 
and 77 = 0.5. It is shown that the threshold value of n T is strongly dependent on the 
value of f2. 
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From Eq.(15), we can find the threshold values of the parameters Q, fix, T and 77, 
beyond which there is not entanglement in the steady state. Some simply inequalities 
can be derived as follows 

°< fi <^ = T + V + nT T ' (16) 

< n T < n Tc = ^ - 1. (17) 

The gray area in Fig. 6 depicts the region where the steady state of the two atoms is 
entangled in the case with T = 0.1 and 77 = 0.5. In Fig. 7, the concurrence C st is plotted 
as a function of the intensity of the thermal field with V = 0.01 and rj = 0.5 for 
four different values of g 2 /A. In Fig. 8(b), we show how two atoms initially in various 
different product states would eventually evolve into the entangled steady state in the 
presence of the external noise driving one of the atoms. Otherwise, in the absence 
of the intense enough external noise, two atoms firstly become entangled due to the 
dipole coupling induced by the virtual photon process, then, they rapidly lose the en- 
tanglement, as shown in Fig. 8 (a). 

Finally, we attempt to discuss the nonlocality of two atoms in the steady state. 
The nonlocal property of two atoms can be characterized by the maximal violation of 
Bell inequality. The most commonly discussed Bell inequality is the CHSH inequality 
[13,14]. The CHSH operator reads 

B = a ■ a <g> (b + V) ■ a + a 1 ■ a <g> (b - 9) ■ a, (18) 

where a, a', b, b' are unit vectors. In the above notation, the Bell inequality reads 

< 2. (19) 

The maximal amount of Bell violation of a state p is given by [15] 

B = 2yJ\ + X, (20) 

where A and A are the two largest eigenvalues of T p T^. The matrix T p is determined 
completely by the correlation functions being a 3 x 3 matrix whose elements are 
(T p ) nm = Tr(pcr n ® a, m ). Here, <J\ = a x , (X 2 = &y, and 03 = a z denote the usual 
Pauli matrices. We call the quantity B the maximal violation measure, which indicates 
the Bell violation when B > 2 and the maximal violation when B = 2^/2. For the 
density operator p st in Eq.(13) characterizing the steady state of two atoms, A + A can 
be written as follows 

A + A = 4| P y 2 + max[4| P y 2 , (p s u + p s u - p s 22 - p* 3 ) 2 ]. (21) 

Recently, Verstraete et al. investigated the relations between the violation of the CHSH 
inequality and the concurrence for systems of two qubits [16]. They showed that the 
maximal value of B for given concurrence C is 2^/1 + C 2 , which can be achieved by the 
pure states and some Bell diagonal states. If the given concurrence C is larger than 



7 



the minimal value of B is 2y/2C, which can be achieved by the maximal entangled 

mixed state. Furthermore, the entangled two qubits state with the concurrence C < ^ 
may not violate any CHSH inequality, even after all possible local filtering operations, 
except their Bell diagonal normal form does violate the CHSH inequalities [16]. So, 
it is not difficult to understand the following results. Our calculations show that not 
any violation of CHSH inequality will be found in the steady state even though the 
steady state is entangled. Moreover, the stochastic-resonance-like behavior can not be 
observed in the Bell violation of two atoms during the evolution, and the stronger the 
noise intensity, the more rapid the Bell violation disappears, which is shown in Fig. 9. 



V. CONCLUSSION 

In this paper, we investigate the problem of generating entanglement when only inco- 
herent sources are available and show that controllable entanglement can arise in this 
situation. We show that, if two atoms are simultaneous driven by two independent 
white noise field with the same intensity, the entanglement between them is suppressed 
and eventually completely destroyed by the noise. However, in another case, in which 
only one atom is exposed in white noise field, the steady state entanglement of the 
two atoms is non-monotonic function of both the intensity of noise driving field and 
the spontaneous decay rate. A double resonance behavior emerges. Moreover, the 
threshold value of the spontaneous decay rate, below which there is not any steady 
state entanglement, increases with the intensity of noise field. 
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Figure Caption 



FIG.l the concurrence C s (t) is depicted as a function of the time t and the intensity 
of the thermal field ny with g 2 / A = 0.2 and T = 0.01. 

FIG. 2 the concurrence C s (t) is plotted as the function of the time t and the coupling 
constant T of the atoms and the external field with g 2 / A = 0.2 and two differ- 
ent values of effective photon number tit of the thermal field, (a) tit = 0; (b) 
fix = 0.3. 

FIG. 3 the concurrence C st is plotted as a function of the spontaneous emission rate 
rj and the intensity of the thermal field tit with g 2 / A = 0.2 and T = 0.1. 

FIG. 4 The concurrence C st depicted as a function of the spontaneous emission rate 77 
and the coupling constant T with g 2 / A = 0.2 and n T = 2. 

FIG. 5 The concurrence C st is plotted as a function of £7 and the intensity of the ther- 
mal field n T with T = 0.1 and rj = 0.5. 

FIG. 6 This figure depicts the region where the steady state of the two atoms is en- 
tangled in the case with T = 0.1 and 77 = 0.5. 

FIG. 7 The concurrence C st is plotted as a function of the intensity of the thermal 
field ut with V = 0.01 and rj = 0.5 for four different values of g 2 / A (from top to 
bottom, g 2 /A = 0.49, g 2 / A = 0.50, g 2 / A = 0.505 and g 2 / A = 0.51). 

FIG. 8 The concurrence C is plotted as a function of the time t with T = 0.1, 
g 2 1 A = 0.2 and i] = 0.5 for two different values tit = 10~ 6 (a) and tit = 2 
(b) of intensity of the thermal field, and for three different initial states: (Solid 
line) |0)i <g) |0> 2 ; (Dot line) |l)i <g> |0> 2 ; (Dash Dot line) |0)i ® |1) 2 (The values of 
the solid line in (a) is too small to be seen). 

FIG. 9 The maximal violation B is plotted as a function of the time t with T = 0.01, 
g 2 1 A = 0.2 and rj = 0.01 for three different values of %, nx = (Solid line), 
ht = 0.5 (Dot line), nx = 1 (Dash Dot line). Two atoms are initially prepared 
in |l)i ® |0) 2 . 
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(a) 




Figure 2: the concurrence C s (t) is plotted as the function of the time t and the coupling 
constant T of the atoms and the external field with g 2 /A = 0.2 and two different values 
of effective photon number tit of the thermal field, (a) tit = 0; (b) tit = 0.3. 
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Figure 3: the concurrence C s t is plotted as a function of the spontaneous emission rate 
T) and the intensity of the thermal field nx with g 2 /A = 0.2 and T = 0.1. 




Figure 4: The concurrence C s t depicted as a function of the spontaneous emission rate 
i] and the coupling constant T with g 2 /A = 0.2 and nx = 2. 
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Figure 5: The concurrence C st is plotted as a function of Vt and the intensity of the 
thermal field ut with T = 0.1 and rj = 0.5. 
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Figure 8: The concurrence C is plotted as a function of the time t with T = 0.1, 
g 2 /A = 0.2 and 77 = 0.5 for two different values tit = 10~ 6 (a) and tit = 2 (b) of 
intensity of the thermal field, and for three different initial states: (Solid line) |0)i<8>|0)2; 
(Dot line) |l)i <g> |0) 2 ; (Dash Dot line) |0)i ® |1) 2 (The values of the solid line in (a) is 
too small to be seen). 



17 



B 




5 10 15 20 

t 



Figure 9: The maximal violation B is plotted as a function of the time t with T = 0.01, 
g 2 /A = 0.2 and r\ = 0.01 for three different values of ut, n T = (Solid line), tit = 0.5 
(Dot line), ny = 1 (Dash Dot line). Two atoms are initially prepared in |l)i <g> |0) 2 . 
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